Conductivity in glass phases of disordered granular superconductors in magnetic fields 



Ryusuke Ikeda 
Department of Physics, Kyoto University, Kyoto 606-8502 
(Dated: June 8, 2009) 

The electric conductivities in glass phases of three-dimensional (3D) granular superconductors in 
magnetic fields are examined based on a disordered quantum Josephson-junction array. Contrary to 
a recent argument of a glass phase with metallic response, a correct inclusion of an Ohmic dissipative 
dynamics always leads to a glass phase with divergent dc conductivity. With no dissipative term, 
a metallic glass phase is obtained irrespective of the range of correlation of quenched disorder, i.e., 
even in the so-called Bose-glass phase with experimentally vanishing resistivity. 
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T. INTRODUCTION 

o i 

O I Recently, it was argued based on a calculation of the conductivity of two-dimensional (2D) disordered Josephson- 
!h ' junction arrays [E0] that a glass phase peculiar to disordered granular superconductors, called as a phase glass (PG), 
Or is not a superconducting but metallic phase with a nonvanishing resistance at low enough T. It is quite important to 
judge this surprising argument correctly in order to clarify whether the resistivity data suggestive of the presence of 
t • | an intermediate quantum "bose metal" phase 0,0] are intrinsic or not. Subsequently, the argument [3,13] of the glass 
phase with a nonvanishing resistance was straightforwardly extended to the 3D case and the systems under nonzero 
magnetic fields, and thus, the so-called vortex-glass (VG) phase [H resulting from a point-like quenched disorder at 
nonzero temperatures was argued Q to be also a metal phase, contrary to the conventional wisdom. [|| However, 
their argument seems to be based essentially on their assumption of nondissipative bare dynamics of the phase field 
0. The neglect of Ohmic dissipation is not acceptable in discussing, at least, the VG phase at nonzero temperatures. 

In the present paper, we show that inclusion of an Ohmic dissipative dynamic term always leads to the conventional 
picture that any glass phase of granular superconductors in nonzero fields is a superconducting phase with vanishing 
| dc resistivity in the plane perpendicular to a uniform magnetic field H || z. In addition, we point out that, with no 
J> • dissipative dynamics of 9, even the so-called Bose glass (BG) phase created by a line-like (columnar) disorder parallel 
to z becomes a metal phase. Based on the experimentally well-accepted fact that the BG phase is superconducting, 
we argue that the argument 0, 0] of the metallic (nonsuperconducting) VG phase based on the use [H, 0, 0] of 
nondissipative dynamics has no foundation acceptable physically. These results are consistent with those of the 
analysis [|[ from the vortex liquid regime based on the Gaussian glass fluctuation. In contrast, the authors of 
Ref.P, 0] have not studied consistently the regime in nonzero temperatures just above a glass transition. 

o' 

03 ' II. CONDUCTIVITY IN GLASS PHASE DUE TO UNCORRELATED DISORDER 



We start from the hamiltonian 



W « = a T,(-^\ - E ^ C ° S (^ - Or + ,a), (1) 



describing a Josephson junction array with a charging energy 2a on each grain, where r denotes the coordinate of 
each site (i.e., grain), /i is the unit vector pointing to possible nearest-neighbor directions, a the lattice constant of 
the simple cubic or square lattice, and 9 r is a phase operator on the grain at r. The model may be extended to 
a more general one including effects of possible dissipation on each grain and of electromagnetic fields. The most 
straightforward method of performing this is to express the model ([T]) into the corresponding quantum action 



S=- dr£ 



exp[i(6» r (r) - r+M (r) - e*A c ^(r) - e*6A ll {r, t)) ] + c.c 
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in the unit % = c = 1, where (3 = 1/T, e* is the Cooper-pair charge, A ext ^ denotes the line-integral of an external 
gauge field over a single bond in the /z-direction, and 8A^{t) is the corresponding gauge disturbance introduced for 
obtaining the conductivity in the /^-direction. Further, 
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is the action corresponding to the charging energy, i.e., the first term of eq.(l). 
Note that the dissipative (last) term of eq.@ is expressed as 



Sdi^/r'EE^MI^HI 2 , (4) 

r uj 

where $ r (^) is the Fourier transform of 

$ r (r) = exp(i0 r (r)). (5) 

That is, eq.Q is nothing but the familiar dissipative term, written in the phase-only approximation, in the time- 
dependent Ginzburg-Landau model. 

Hereafter, let us proceed to rewriting the action S into a form convenient for a field-theoretical method. Further, a 
point-like (uncorrelated) disorder will be assumed in this section. A quenched disorder in the system is incorporated 
into a randomness of = J* with a nonzero real mean Jo, i.e., J M = Jq > 0, and a Gaussian distribution 

(J^^i) — ./o)(^-^( r 2) — Jo) = 4£ ri . r2 J 2 /d, where d is the space dimension. These relations may be regarded as being 
due to a random gauge field a p defined by J p (r) — J oc exp[ia M (r)]. The free energy F = — /3 -1 lnZ will be expressed 
in terms of the replica trick as F = —f3~ 1 (Z n -l)/nin?i-> +0 limit. The averaged replicated partition function Z n 
is given by 

= Zq < exp(-S/ - S g ) >o, (6) 
where Z is the partition function of Sq, < >q denotes the ensemble average on Ei< a <n^(^°'): an< ^ 

S ^ = -EE / dr^cos(e*SA,(r) + e*A cx ^r)-e^(r)+ei% a (r)) 

a=l r,/i " , ° 

- i fdr f dr'^^J 2 cos(e*[<5^(r)-M M (r')] 

J J a,b T,n 

+ 9^(r) - 0<V) - ei% a {r) + 6i% a (r') ). (7) 
Before proceeding further, 5/ will be rewritten in the form 0] 
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5/ = const. -dJo^^E^H^-^H + ^E^W-W)*'^) ( 8 ) 
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for the cubic or square lattice in d-dimension. Equation (9) is valid up to the lowest order in v\uj\/Jq and in the 
gauge-invariant gradient on the lattice [9( accompanied by the gauge field A cx + <5A(r). Then, by introducing the 
conventional SC order parameter ip( a )(r, r) and the PG order parameter ql ab ^ (t\, t-x) = (<Zr (t2,Ti.))*, Z™ becomes 



7n r 

= = / PV (a) ^(V- (a) )*^ (Qb) exp(-5e ff ), 
z J 
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where 
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and Z) M denotes the gauge-invariant gradient on the lattice accompanied by the gauge field SA(ti) — <5A(t 2 ). Performing 
the cumulant expansion in powers of and ip^ in the logarithmic term and replacing (lpj q ab (ri, r 2 ) by Q afc (ri, t 2 ) — 
C<5 aj f,(5(Ti — t 2 ) with a constant C, we finally obtain the following effective Landau action 



tS eS ^,Q;SA) = /5[/^E 
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Further, the r-sum was transformed into the space integral, and instead, a short length cut-off a was introduced. Note 
that the 5A-dependent term arises from a term bilinear in g( ab )(Ti, r 2 ), and a Q( aa )-linear term carries no SA. 

We note that, although the dissipative term in eq.JT]) directly appears only in the term quadratic in ip of the effective 
action, it also affects the dynamics of the glass field through the coupling (uty-) term between the two fields. 

The above expression of the action is of the same form as those in other works [J 0] . As explicitly examined in 
Ref.Q, the coefficients r, u, t, k, c^, d^, and are positive, while 

1 T T rn T//\2 

?>,o 
t 



_1__ Jo JoT_ Jo. . _, _! 
Ad 2a 2a 2 4a\\2J 



(16) 



~i 2 Jo 
a = a — — . 

Ada 

Below, the Fourier transform of the glass field Q( ah \r) is defined, by following Ref . [l(| ■ as 

Q^ ah \v-T llT2 ) = ^) + ij5 w e- Mn ^^ 



^ 2 E SQffltWt 
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where 

=q(l-S a ,b) +qS a ,b (18) 

in the replica-symmetric approximation adopted in the previous [2, 0, E3] papers. 

Now, we examine the conductivity in PG at low enough T by, as in other works [ESS E3|, treating the PG order 
parameter field in the mean field (MF) approximation where SQi ab )(r) — 0. We note that, in eq. (fT5)) . the ^-field (i.e., 
SC fluctuation) couples to the glass field Q l - ab ) in a bilinear form (ip( a ) )* i/jW of which the diagonal (a — b) component 
is nothing but the SC fluctuation contribution to the entropy density irrespective of the presence or absence of the 
applied field. Thus, a result on the conductivity in zero fields (H — 0) is trivially extended to the H > case. In other 
words, the metallic conductivity [3, 0] in PG, if correct, would also affect the H-T phase diagram at low enough T. 
Bearing this in mind, we will discuss, for the moment, results in the simpler H — case with A cx = and concentrate 
on results valid for both the H = and H ^ cases. Our analysis is different from that in other works [l], Q in 
that we take account of a dissipative dynamics correctly and of the presence of the SC fluctuation consistently in 
determining the MF solution of Q( ab \ As shown below, the coupling between the SC fluctuation and the MF Q( ab ) 
makes the conductivity in the PG ordered state not finite but divergent and hence, makes the PG a superconducting 
phase. 

Deep in the PG, it is sufficient to keep the ^-fluctuation in the Gaussian approximation. In H ^ case, the 
corresponding approximation may be valid above H C 2(0) and at low enough temperatures Below, the replica 

symmetry will be assumed even for the static (u) = 0) components on the basis of the argument in the literature [l(| 
that taking the replica symmetry is justified in low T limit. A breaking of replica symmetry cannot qualitatively 
affect our main conclusions given below, because it would be accompanied by an independent parameter such as a 
coefficient of a quartic term on Q ab in the effective action. Then, using eq. lfTT)) and neglecting 5Q( ab \ the variational 
equation = lim„_^ + o"- 1 d 'Z n / 'dQ( ab ) (t\, rg) takes the form of the following three equations: 

/e -1 (w a + r) - K Dl+u(q + /T 1 ^^) (19) 

for nonzero u>, 

K- l r - np 2 {q 2 -q 2 )+u(q + /T^A,) (20) 

- ^^cf (k;0) = o, 



2 K (3 2 q(q -q)+ m / G< J (k; 0) = 0, (21) 

where J k denotes J d d (ka) / (2Tr) d . These are obtained as the variational equations on D u , q, and q, respectively. Here, 
we have expressed the SC fluctuation propagator (k;a>) = /3 _1 < {ipf" (k; u))*t/)( (k; u>) > in d-dimension in the 
form G$(k;a;) = 6 a ,b(l - <W)G^(k; u) + <L, [<5 a , b G^(k; 0) + (1 - *„ )6 )G$(k; 0) ], where 

G ^ w ) = ta i n — br^ — w (22) 

f> + dip\u)\ + Ctyor + ta z k z — w^D^ 
G^(k;0) = ^( 9k (Aq)) 2 , 

G^(k;0) = 5k (A g )+G^(k;0), 
5k(Aq) ' 



r^j, + td 2 k 2 + w^pAq' 
and 



Aq = 0{q-q). 



(23) 
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Noting that, when q > 0, eq. (|21j) becomes 



we easily find that the only physically meaningful solution of the PG order parameter is given together with eq. ([24 
by 



D u = -Aq-K- X p^\u;\y 2 , (25) 
q = -{3- 1 ^T~D U +u- 1 { K {Aq) 2 J <? k (Ag) 



- « -1 r)> 

where 

_( K^w^^ 1 / k (g k ) 2 \ 1/2 

The above form of D u is valid up to Od^l 1 / 2 ). When the dissipative term in G^(k, u>) is absent, the — l^l 1 / 2 term 
in eq.(24) is replaced by a — |w| term. In — > limit where ip and Q fields are decoupled, the above MF solution 
reduces to the one — with q > and ip = if higher order terms, omitted in eq. ([!?!)]) . are kept. Note that 

the |w| 1//2 -term, arising from the dissipative dynamics in eq.©, of D u was brought by the ^-fluctuation. Consistently, 
this term appears in G^(k;u) and makes the dynamics of ^-fluctuation in the PG state sub-Ohmic. 

The analysis in zero field case given above is applied to granular superconductors in nonzero fields in the following 
way: As in eq. (pH|) , effects of SC fluctuation on the glass order are included merely in a form of an integral over the 

momentum k of (<?k)™ 1 or G^(k;w). In nonzero fields (H > 0), the -0-field is decomposed into the Landau levels, 
and one has only to replace the kj_-integral and k 2 ^ by a summation over the Landau level index I and \e*\H(2l + 1), 
respectively, where e* is the charge of Cooper pairs. Thus, it is easily understood that eq. (f25)) and the divergent 
conductivity given below are also valid in H ^ case if such replacements are performed. 

Further, we note that the dissipative term, eq.(4), may be alternatively incorporated as follows: To simplify our 
description, any spatial variation of the glass order parameter q will not be considered so that the gradient may 
be negligible. Further, the dissipative term will be treated here in the form not of Sf but of S g = Sd ls + S g . Explicitly, 



■- \ j j „ h 



(r,r') 



ln(T T exp 



^EEMI^HP + ^E [drf dr¥; b \r,r')q^(r',T) + c.e. 

ui a a,b J J 



(27) 



where ^^ b \r, t') — $£°^(T)($r^(T'))*. Here, let us focus on the replica diagonal terms. Then, expressing </ aa ' (r, r') 
as exp(— iw(r — r')) and replacing d^ — by d u , the above expression is rewritten as 



ln(r r exp 



\Y,JdrJdr% {ab \r,r') q ^\r',- 

a,b 



(28) 



Except the first term, the above expression is the same form as that in the nondissipative case. Contrary to the claim 
0, the first term is not lost after the replacement q^ (r, r') — > Q^ ah \r, t') — C<5 Q ,b(5(r — r'). That is, the \lu\ term is 
added to the first (lo 2 -) term in eq.(18). Then, it will be obvious that the variational solution of Q( aa *> takes the form 
oc — l^l 1 / 2 as its leading uj dependence at low |u;|, just as in eq.(24) (see the sentences following eq.(25)). In Ref.[3], 
the first term was erroneously omitted, and consequently, the argument of the metallic glass phase was not changed. 

Now, let us examine the conductivity a following from the solution, eq.(25), in terms of Kubo formula. When the 
glass fluctuation SQ^ is neglected, the conductivity following from the present model consists of two parts. One is 
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the direct contribution from the SC fluctuation and, as usual, has the form 



4(e*) 2 a 2 



(29) 



Wi a, b 



[Gil ) (k;c 1 )-G^(k;^i+^)] 



iw->f2+iO 



where the n — > limit is taken at the end. In the disordered (i.e., normal) phase in which q — 0, this expression was 
examined previously and shown to vanish in T — > limit (l3j although the bare dynamics is dissipative. In the PG 
phase, by substituting eqs. (|22)) with (|25|) into eq. ([29)l . er^ becomes 



Recr v ,(fi ->■ 0) 



2V2q(e*) 2 p l /) ^ 
/s|n|Va 



fc 2 (5k) 3 



(30) 



for d = 2, implying a divergent contribution to the dc conductivity arising from the sub-Ohmic dynamics of the SC 
fluctuation in the PG phase. Note that the T — > limit was not taken in obtaining eq. (|30p . 

To corroborate that the total conductivity is also divergent, let us examine another contribution 



cr PG (iw) 



— — I — I Urn - 7 / drie 
\u\d\aj ^on^/p 



iw(ri— T 3 ) 



(31) 



dr 2 S( n - r 3 )\Q^(n,T 2 )\ 2 - \Q (ab, (n - r 3 )\ 



(aft). 



to the conductivity which arises from the terms (see eq. fTl")) ) quadratic in Q( ah ^ and dependent on SA. By substituting 
eq. (fT7)) into eq. (j3Tj) , it becomes 



a PG (iLo) = 



1 /e* 
M d \ a 



(1 - <L,o) 







E q (ab) q {ba) + /r 1 E D "* ( D «i - D ^+») 



(32) 



a. . b 



l3- 1 (2D OJ -Do)D -2q^ £>, 
1 



(i - « w ,o) - E 9 (Qb) 9 (fca) + V £^ - ^^+-) 2 + 2 <^w - 



where the replica-symmetric expression (|18l) of g( ab ) was used in obtaining the expression following the second equality. 
Consequently, we have q^ ab ^ = q, and the sum q( ab )(f ba ) term is absent in n — > limit. Then, let us turn to the terms 
including D u . The real part of the term under the wi-summation is positive and nondivergent after the substitution 
ic<j — > £1. It can be easily seen in terms of the identity 



dx- 



(33) 



In contrast, according to eq. ([23)) . the term proportional to q is positive and, as well as eq. ([3"0"]) . becomes divergent like 
|f2| -1 / 2 in dc limit. The above results altogether imply that the PG phase is superconducting. We stress that this 
result has been obtained by including the coupling between the SC (ip-) fluctuation and the glass order parameter 
Q(ab) in obtaining the MF -solution. 



III. EXTENSION TO GLASS PHASE DUE TO LINE-LIKE CORRELATED DISORDER 



Here, our analysis of the conductivity will be extended to the case of the so-called BG state which occurs due to 
line defects persistent along the applied field || z. Then, we assume the starting model eq.(2) to, in turn, have the 
correlation J M (r) J*, (r') = (4J 2 /d)^ ljM '<5 x ,x' of the random Joscphson coupling, where r = (x, z), r' = (x', z'), and 
x (x') is the in-plane component of the coordinate r (r'). To point out our main point that, with no dissipative 
dynamics, the physically acceptable response properties in a glass phase do not follow, we only have to focus here on 
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the simple case with no coupling to the -0-fluctuation and with no dissipative phase dynamics. If incorporating the 
gauge field 8 A _L z necessary in examining the conductivity perpendicular to the field, the replicated action related 
to the glass field or Q( ab *> takes the form 



a.b x z,z' 

x exp(i((L4 M (r) - %(r'))) + c.c. 

-J 2 f , f , /\-^ \ ^ l^-(ba), , 



$ { * a \z', z; t', t) $i ah) (z + na, z' + fia; r, r') 



dr / dr'Y: E 

a.b x,z,z' 



D„ ■ D 



2d 



C(v';r,r') + cc. 



(34) 



where $i &a) (z' , z; t', t) = ^ {z 1 , r') } (z, t))*, and $i a) (z,r) = cxp(i6»( a ) (x, z; r)). Then, by assuming the cases 
with the glass field independent of x, the effective action corresponding to eq.(ll) and with no V'-field becomes 

J~ 



x ' a, b z.z ! 

- E ln(T T exp[iEE j dT j dr'^ b \z,z':r,r')(l-^\^(z',z;r',r) + 

x L a, b z,z' 

= E[-EE / dr%^G(r-r')fl-^) q ^(z',z;r',r) 

a, 6 j=l Zj J ^ ' 

1 6 /" 

3 E II E / dT 3 G ( T i - t 6)G(t 2 - t 3 )G(t 4 - T 5 )8z 1 ,z 6 5z 2 , Z3 Sz i ,zAf ) 934° 956° + ' 

a,b,cj = l Zj 

where qmn = q( ab \z m , z n ; r m , r„). After the transformation into the Q( ab ) -field, this action becomes 

d 2 



(35) 



tSeS = E EE K 1 f dTl f dT 2 S ( T l - T 2) S z 1 ,z 2 (g^ 



l9t2 4d dz\dz2 



■-zi,Z2 a 
3 



Q&°(x) 



a,6,c j — 1 ' 



+ 



E E j**J E (H(n)-H(^)) 2 QlfC 



a, 6 ^i}^2 ' 



Q 



(36) 



This form of the effective action implies that the nonlocality in the z-dircction || H appears in parallel with that in 
the imaginary time. Hence, if the glass fluctuation is negligible, the glass field should take the form 



Q {ab) = qW{ ZuZ2 . TuT2 ) = q W +p- 1 L- 1 Y,DUk z )S a . b cM^z(z 1 -z 2 )-iLo(r 1 -T 2 )) 



(37) 



iv, k 



as an extension of that in the point disorder case. Then, it is straightforward to verify that, in the replica symmetric 
approximation, 



lab = q, 



(38) 



This form of the Q/^-'-solution indicates that, as far as focusing on the Kubo formula of the conductivity, contribu- 
tions to the conductivity from the k z = components of Q^ aa '> should be the same, in the BG phase, as those in the 
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2D VG point phase at T = 0. To see this, let us first focus on the counterpart of eq.(32). In the case of line disorder, 
it is easily seen that, using eq.(37), the counterpart of eq.(32) becomes 



JBG) 

7 PG 



1 



lu\ d 





iim I f dT ie iu(Ti - T ^ 



(39) 



Z1,Z2 a,b 

dr 2 6(n-T 3 )\Q^\ 2 



\Q 



(afc),2 
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and that the terms proportional to q (see the final term of eq.(33)) in eq.(40), which are the only contributions there 
leading to a nonzero conductivity in T — > limit, consist only of the k z = contribution. Thus, just as in the 
nondissipative point disorder case, a nonzero positive contribution to the conductivity follows from eq . (f3"9"| . Further, 
the contribution, eq.(31), related to the SC fluctuation also gives a positive contribution to the conductivity stemming 
from D UJ (k z = 0) and is expressed by 



JBG) 



(iu)) 



d^2q(t 



a =o) 



(40) 



Again, this SC fluctuation contribution is also finite in the nondissipative case. Therefore, the results on the conduc- 
tivity in BG phase are qualitatively the same as in VG phase. By repeating the analysis mentioned above in the case 
with dissipative dynamics, it is easy to verify that the contributions to the conductivity corresponding to eqs.(40) and 
(41) are divergent in the dc limit |fi| — > in the case with dissipative dynamics. 



IV. SUMMARY 

First, we have shown by applying a theory of quantum spin-glass to the quantum granular superconductor, 
modelled by eqs.{T|) and (JSJ, that, even if assuming the replica symmetry, a PG ordered state peculiar to granular 
superconductors is characterized by zero resistance when the Ohmic dissipative phase dynamics is incorporated. This 
is consistent not with the argument [l[ of a metallic response in this phase but with the conventional picture [TH [l2j 
of a superconducting glass existing in H > H C 2(0) (see Fig.l). Elsewhere, we have shown [8] that, even without the 
PG order, the fluctuation of the PG order parameter plays the role of pinning disorder inducing a VG instability at 
T > 0, at which the resistivity continuously vanishes, in 3D and nonzero field case. The present result is consistent 
with results of the approach [8| from higher temperatures. 

We have also examined the case with the glass (BG) phase due to line-like disorder in addition to the familiar VG 
case due to point disorder in order to corroborate that the argument [7J that the VG is a metal is false even from 
the viewpoint favoring consistencies with experimental facts. It was found that, even in the case with line defects, 
the conductivity in the resulting BG phase also becomes nondivergent when the phase dynamics is nondissipative, 
while it is divergent in dc limit when the Ohmic dissipation is introduced. Thus, the argument in Ref.0] based 
on the nondissipative model that the VG is a metal is equivalent to concluding that even the BG is also a metal, 
contrary to the well accepted experimental fact that the BG has zero dc resistance, and hence, has no reliable basis 
both theoretically and experimentally. The result that, in both cases of point-like and line-like disorders, the glass 
phase becomes superconducting in the presence of dissipative dynamics indicates that there is no theoretical basis 
of resulting in an intermediate bose- metal phase in disordered granular systems at low T limit. When trying to 
understand superconducting resistive behaviors upon cooling, inclusion of a dissipative dynamics is essential even in 
the quantum regime 
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FIG. 1: Schematic mean-field H-T phase diagram of a disordered granular superconductor. The dashed curve implies 
microscopic H c2 (T) curve for a single grain. 



